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I. INTRODUCTION 



Semiclassical methods have proved to be very useful in the investigation of systems with 
many degrees of freedom, especially in atomic and molecular dynamics^"^. Moreover, the 
semiclassical approximation has also been an important theoretical tool in studying the 
connection between the classical and quantum theories, particularly in fundamental topics 
such as chaos and open quantum systems'"'. 

The semiclassical propagator in the coordinate representation was first derived by Van 
Vleck'^ at the beginning of the last century. However, this fundamental result has two re- 
markable characteristics that considerably hinder its practical application. First, the Van 
Vleck propagator is determined by classical trajectories subject to boundary conditions. In 
general, the search for these specific solutions is quite complicated, particularly in multidi- 
mensional and chaotic systems. The second major problem is the appearance of focal points, 
which are responsible for divergences in the semiclassical approximation. 

A different line of research, concerned with the difficulties caused by focal points, led to 
the development of semiclassical propagators in the representation of the harmonic-oscillator 
coherent states""^-. Although it has been found that the focal points still persisted, this 
alternative approach has demonstrated some evident advantages over the coordinate and 
momentum representations, including an immediate visualization of the system over the 
full phase space. Nevertheless, new problems have emerged, such as the duplication of 
the phase space, resulting from the apparent overdetermination of the classical equations 
of motion. Furthermore, not all classical trajectories in the extended phase space, while 
correctly satisfying the boundary conditions, correspond to semiclassical propagators with 
physical meaning ' . Therefore, it is necessary to establish effective rules for selecting the 
proper contributions to the semiclassical dynamics. 

In the last decades, many different techniques have been proposed in order to solve the 
recurrent problems in semiclassical propagation Most of these methods are based on 
the concept of initial value representation, in which the system dynamics is determined only 
by initial conditions, avoiding the search for boundary-valued trajectories. 

Recently, Aguiar et al. presented a new approach to the semiclassical propagator of 
the harmonic-oscillator coherent states, which combines the unique resources offered by 
the trajectories in a doubled phase space with the plain advantages of an initial value 
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representation'^^. Moreover, they demonstrated that practical and simple rules for selecting 
the contributing trajectories can produce very accurate results. 

The procedures developed in the present paper are similar to those of Aguiar et al., but 
generalized to a subclass of the SU(?2) coherent states. These states constitute an ideal 
setting to study the bosonic dynamics for a fixed total number of particles in n modes. In 
this paper we also propose a new prescription for the selection of contributing trajectories, 
which we designate as a heuristic filter. For a detailed derivation of the SU(n) semiclassical 
propagator we refer the reader to a recent work of the present authors'^'. 

The remainder of the paper is organized as follows: in section II we develop the semiclassi- 
cal propagation method based on an initial value representation. We start with a brief review 
of the SU(n) coherent states, in which we introduce fundamental aspects of the adopted no- 
tation. Then, we present the SU(n) semiclassical propagator, followed by other important 
definitions, such as the effective classical Hamiltonian, the classical equations of motion and 
the doubled phase space. Next we reformulate the semiclassical approximation in terms of 
a set of initial conditions and a heuristic filter of trajectories. At the end of the section, we 
describe the procedure used for calculating semiclassical mean values of observables, based 
on the phase space representation of states. Section III presents an application of the SU(2) 
and SU(3) semiclassical propagators. As an example, we consider a simplified model for 
the dynamics of a Bose-Einstein condensate in a triple-well potential. In this context, we 
introduce the classical approximation, which provides a reference for comparison with the 
semiclassical results. Also, we discuss the accuracy of the semiclassical propagation in non- 
linear and predominantly linear dynamical regimes, by contrasting the approximations with 
exact quantum calculations. Finally, in section IV we present our concluding remarks. 

II. SEMICLASSICAL PROPAGATION METHOD FOR SU(n) 

A. SU(n) coherent states 

The coherent state related to the fully symmetric irreducible representation of SU(n) for 
identical bosons is given by^^: 



w) 



mi+m2+...+mn=N 



E 
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mi\m2\ . . . mn\ 



N\ 
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mi,m2, . . . ,mn) 
(1 + w*w)^ 



(1) 
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where {|mi, m2, . . . , m„)} is the usual basis of the bosonic Fock space for n modes 
and particles, such that nij is the occupation in the j-th mode. The vector w = 
{wi,W2, ■ ■ . ,Wn-i)^, with (n — 1) complex entries, parametrizes the entire set of coherent 
states. 

Although normalized, the coherent states in (1) are not orthogonal'^^: 

{w \w) = ^ -w- (2) 

(1 + W *W ) 2(1+ W*W) 2 

However, due to the overcompleteness of the coherent states, we can write the following 
diagonal resolution for the identity in B^: 

n-l 



/a(n) dim(B^) -p-i- , 
dLi(w ,w) \w){w\ = 1: da(w ,w) = — r — I I a Wj\ 



(3) 



where d'^Wj = dxjdyj, with Xj = Re (wj) and yj = lm{wj). Note that the normalization 
factor in (3) can be divided into a{n) = ^"^n-l' , which is independent of the total boson 
number, and dim(B^) = ^^^^^I^y' , the dimension of the accessible Hilbert space. 

B. SU(n) semiclassical propagator 

The quantum propagator in the SU(n) coherent state representation is defined as the 
transition probability between the initial coherent state \wi) and final coherent state \wf) 
after a time interval r "': 

K{w},wf,T) = {wfle-'^'^lw,). (4) 

After recasting the above propagator as a path integral, we can perform its semiclassical 
approximation, which consists in expanding the action functional to second order around a 
classical trajectory. The result of this derivation is given by^^: 



1 + w(t)w(t) 



2 

det 



dw{0) 



(5) 



_l+w(0)w(0)J [dw{T)_' 
All elements of this semiclassical formula are calculated on a classical trajectory, which 
is solution of the equations of motion'" 

w = - — (1 + WW) [I +w ^w\ -— = -tt,^— 
N ow aw 

^ ^ _ ^ - -OH ^ ' 

w = —il + WW) [1 + w ^w\ —— = it.^:— 
N ow ow 



with boundary conditions 



w{0) 
w(t) 



Wi 



(7) 



w 



In equation (6), H is the effective classical Hamiltonian: 

{w*\H\w) 



T-L(w, w) 



{w*\w) 



(8) 



If the classical equations of motion have more than one solution subject to the same 
boundary conditions Wi and w*j with fixed time interval r, then the correct semiclassical 
propagator between these points is given by the sum of the propagators (5) for each possible 
trajectory. 

Note that the complex vector variables w and w are completely independent, i.e. in 
general w{t) ^ w*{t). This doubled phase space is a direct consequence of the introduction 
of boundary conditions to the equations of motion. If w{t) were equal to w*(t), the two vector 
differential equations in (6) would be redundant and the boundary conditions w{0) = Wi 
and w*{t) = w*j would make the problem overdetermined. Therefore, the duplication of 
the phase space is required to solve the classical equations of motion in the coherent state 
representation. 

The equations of motion (6) are derived by the extremization of the following action 
functional: 



S{w*j,Wi;T) = J L (^w,w,w,w) dt + T{wj,Wi;T); 



WW — WW 



L (w, w, w 



v,w) 



H{w,w) 



(9) 



T{w*j,Wi;T) 



2 1 + WW 
-t^Ln { [1 + w}w{r)] [I + w{0)wi] } . 



The function F, known as the boundary term, is essential in obtaining the classical 
equation of motion subject to the boundary conditions (7). Another quantity introduced in 
(5) is the correction term to the action^'^: 



Tr 



— (l— 
dw V dw 



dw \ dw J 



dt 



(10) 



where the matrices ^ and ^ are defined in equations (6). 



The last ingredient required in the formula (5) is the tangent matrix M, governing the 
dynamics of small displacements around the classical trajectory, defined in block form by 





Mn(r) Mu{t 
^M2i(r) M22(r 

Notice that 

M22{T) - - 




(11) 



(12) 



dw{0) idw{T) 

and, therefore, the block M22 is the inverse of the matrix whose determinant appears in the 
semiclassical propagator. A focal point in the variables w corresponds to a zero value of 
det M22(t) and, consequently, to a divergence in (5). 

The tangent matrix can be calculated as solution of a system of differential equations 
subjected to initial conditions. Using (6), we obtain 




_-_d_ \cdU] _i_d^ 

dw IS dw J dw [y dw J 

■_d_ \c&H] [£^1 

dw V ^ dw J dw L ^ dw \ 





(13) 



Substituting the definition (11) in (13), we find 

M(t) = M(t)M(t). (14) 

with initial conditions 

M(0) = 1. (15) 

However, note that the matrix ]R(t) is calculated on the classical trajectory, which in 
its turn is subject to boundary conditions. Also notice that the differential equations (14) 
couple the blocks of the tangent matrix exclusively in pairs. Therefore, we need to consider 
only the equations of motion for Mi2(t) and M22(t), with initial conditions Mi2(0) = and 

M22(t) = 1. 



C. Initial value representation 

The classical trajectory is the fundamental quantity for calculating all elements of the 
semiclassical propagator. However, finding the classical solution represents a boundary con- 
dition problem, whose analytical or numerical resolution generally exhibits greater technical 
difficulties or higher computational cost than a similar problem subject to initial conditions. 
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Therefore, the development of semiclassical propagation methods based on initial conditions, 
known as initial value representations, is highly desirable. In this section we develop such a 
method for (5). 

First, we use the resolution of the identity (3) to reconstruct a specific propagator from 
an integral over the entire set of propagators with the same initial coherent state: 

K{w},w^;t) = {wf\e-'^^\w,) 

dfi{w* (r) , w{t) ) {wf\w* (r) ) {w* (r) | e-'^^\wi) 



dii{w* {t),w{t)) {w f\w* (t)) K (w{t) , Wi] t) 

Next we consider w[t) as a function of the initial values of its corresponding trajectory: 

w{t) = w{wi, Wi]T); (17) 

where Wi = w{0). Thus, the integrand in the last line of (16) also becomes a function of 
Wi implicitly in w{t). The change of integration variables introduces the following Jacobian 
determinant: 



ra-l 



l[d%{T)= det lll^ l[d'w,{0) = \detM,2{T)fl[d%{0). (18) 
j=l L ^ '^-l j=l j=l 

We should note that the mapping between Wi and w{t) is not injective, due to the 
existence of focal points. However, the determinant of M22 is zero at these problematic 
values of w{t), so that their contribution to the integral is nuU^"'. 

Finally, considering the semiclassical approximation for the propagators in the integrand 
and substituting the expression (18) in (16), we obtain the first form for the semiclassical 
propagator in the initial value representation: 

r^^vr, * A f ^2- dim(B^) |det M22(r) |' , , 

KZ'iw},w^;r)= j d^w, ' {wf\w*{T))K,,{wiT),wf,T); 

(19) 

n-l 

where d'^Wi = Yl d'^Wj{0). Notice that the integrand of (19) is now proportional to 

3 _ 1 

|det M22(t)| 2 , instead of the inconvenient factor |detM22(T)| ^ in equation (5). Thus 
we avoid the potential divergences of the semiclassical propagator corresponding to focal 
points in the variables w. Also note that all quantities in the integrand of (19) are calculated 



2 n~l n— 1 



on the trajectory with initial conditions w{0) = Wi and u'(O) = Wi. Therefore, by calculating 
the semiclassical propagator Ksc{w{wi,Wi]T),Wi,T) for a grid of initial conditions with Wi 
fixed, we obtain the semiclassical propagator w,; r), at the desired arrival point, 

after an integration in Wi. 

However, our scheme to recast the propagator in terms of initial conditions seems to have 
some disadvantages in relation to the original boundary condition problem. At first glance, 
we replaced the calculation of a single propagator by an infinite number of propagators, 
which are calculated for all possible values of Wi. Even though the latter are subjected to 
initial conditions, the large number of propagators in the integration can make this method 
impracticable. But experience tells us that the trajectories with major contribution to the 
integral (19) are associated with values of Wi close to w*. Therefore, the integral (19) is 
usually calculated for a small grid around w*, considerably reducing the number of classical 
trajectories required in a practical application. 

The second problem in the expression (19) is the need to carry out a new integration for 
each choice of the final coherent state, parametrized hj Wf. However, all dependence on Wf 
in the integrand of (19) comes from the factor {wf\w*{T)). Hence, using the identity (2), we 
can perform a multinomial expansion in the numerator of the coherent state overlap, thus 
extracting Wf from the integration sign: 



KT:(w},w.;r) 



m-i + ...+mn=N 



mil . . .m„! (1 + wjWf) 



n-l 



U=i 



-mi ,...,r7i„ 



(20) 



Hence, in order to calculate the semiclassical propagator for an arbitrary final coherent 
state, we need to perform only dim(B^) integrations whose values are independent of Wf-. 



2^711 ,...,m„ ("U^j ) 



dim(B^)|detM22(r)|' JJ [ 



mi\m2\ . . ■ mj. 



_ 1 

2 

(mi,m2,. . . ,m„|e"*^^| 



Ksc{w{T),Wi; t) 



(21) 



The second equality shows that the integrals Tmi,...,m.„ can be rewritten as semiclassical 
propagators between the initial coherent state and a number state, except by a combinatorial 
factor. 
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D. Heuristic filters 



It is well known that some trajectories in the doubled phase space give unphysical contri- 
butions to the semiclassical propagator Therefore, given a grid of initial conditions 
Wi, only part of the resulting classical trajectories participate in the calculation of the inte- 
grals (21). The appropriate contributions can be collected using the heuristic filter defined 



The classical trajectories that violate this condition at time t are discarded from the 
integration for t > t. Note that the only free parameter in the initial value representation is 
A, whose positive value should be adjusted in order to optimize the semiclassical propagation. 

The idea behind this filter is the following: if we write the semiclassical propagator as 
Ksc = e°'~^^^, with a, (3 E H, then the inequality (22) can be recast in the form ^ < |- 
Therefore, the discarded trajectories are those that lead to an abrupt positive change in the 
real part of Ln {Kgc), thus causing the divergence of the absolute value of the propagator. As 
seen in the equation (5), the time variations in a are directly determined by the imaginary 
part of the corrected action {S + 1). However, unlike previously published methods"''"', the 
proposed heuristic filter also takes into account the factor that contains the determinant of 
the tangent matrix. Clearly, the modulus of this factor also affects the value of a, either 
counteracting abrupt negative changes in Im(S' + /) or contributing to the divergence of the 
semiclassical propagator. The inclusion of this aspect in the heuristic filter is an important 
element in the present work, which greatly improved the results in section 111. 

E. Q representation with SU(n) coherent states 

Using the expressions (20) and (21), we can easily calculate the semiclassical propagator 
at any point w of the classical phase space for fixed initial condition Wi and period of 
propagation r. Thus, we obtain a complete description of the system state, known as the 
Husimi or Q representation . In general, the function Q{w*, w) associated with an arbitrary 
state lip) is defined as: 



by: 




(22) 



Q{w* ,w) = {w\p\w) 



(23) 



\Kiw*,w,;T)\'; 
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where p = is the density operator for a pure state and \w) is given by equation 

(1). In the second hne of (23) we assume that {ip) = e~^^'^\wi). Therefore, using the 
SU(n) semiclassical propagator, we can directly construct the semiclassical representation 

Qsc{w*,w) = \Ksc{w*,Wi,T)\'^. 

With the aid of the expression (3) and assuming = 1 we find that: 



Unhke the exact definition (4), the semiclassical propagators (5) and (20) do not preserve 
the norm of the state during its evolution'^*^. Therefore, for a proper comparison with the 
quantum results at time r, we need to normalize Qsc according to the relation (24). The 
normalization of the quantum and semiclassical representations is implied in the remainder 
of the paper. 

In terms of the exact Q function or of its semiclassical version Qsc, we can readily obtain 
the mean of an arbitrary observable O: 



The function O^, which corresponds to the antinormally ordered symbol of the operator 
O, is defined by 



III. APPLICATION AND DISCUSSION OF THE SEMICLASSICAL 
PROPAGATOR 

A. Bose-Einstein condensate in a triple-well trapping potential 

In order to illustrate the method described in previous sections, we discuss here its appli- 
cation to SU(2) and SU(3) coherent states, considering a simplified model for the dynamics of 
a Bose-Einstein condensate in a triple- well potentiaP*^. Assuming that the three wells of the 
trap are identical and equivalently coupled, the Hamiltonian of the model in a three-mode 
approximation is given by: 




(24) 




(25) 




(26) 



3 




(27) 



10 



where aj (at) is the bosonic annihilation (creation) operator related to the single-particle 
state \uj), which represents the ground state of a harmonic oscillator centered on j-th min- 
imum of the trapping potential, for j = 1, 2, 3. The parameters Q and x correspond to the 
rates of tunneling and collision of trapped bosons, respectively. 

Note that H preserves the total number of particles, so that we can restrict our analysis 
to invariant subspaces with fixed A^, denoted by B^. Therefore, the SU(3) coherent states, 
defined in equation (1) with n = 3, are appropriate to study the model. Substituting (27) 
in (8), we obtain the effective classical Hamiltonian: 

v. ^WiW2 +W2W1 +Wi + Wi +W2 + W2 

1 + WiWi + W2W2 , s 



N 1 + wiWi + W2W2 

-t- w'^Wn + 1 

+ X-, 



(1 + WiWi + W2W2Y' 

Then, employing the general formula (6), we find the classical equations of motion for 
the condensate: 

iWj = n{l + wi + ^2 1 - Wj) + 2x- — _ ^ ^ _ ^ 

' 1 + WiWi + W2W2 ^29) 

^/ N / N WjiWjWj — 1) 

-tWj = n{l +Wi+ W2){1 - Wj) + 2x , ^ ^ _^ 

1 + WIWI + W2W2 

for j = 1,2. Using (28) and (29), we can easily obtain the other dynamical quantities 
relevant to the calculation of the semiclassical propagator, such as the Lagrangian L and 
the matrix M. According to the equations (29), the dynamics of the condensate exhibits 
three classical invariant subspaces, described by the following conditions: 

Wi = W2, Wi = W2 (30a) 
wi = 1, wi = 1 (30b) 
W2 = I, W2 = I (30c) 

For simplicity, we limit our discussion to the case (30a), since the three invariant subspaces 
are dynamically equivalent^'"'. Now, we show that the effective quantum dynamics of the 
condensate under the constraints (30) can be approximated by SU(2) semiclassical prop- 
agators. For this purpose, we rewrite the coherent state (1) in terms of bosonic creation 
operators: 



\w} = , 



n-l 



n 



j=i 

(1 + w*w)2 



N 



|0). (31) 



11 



Then, we apply the condition (30a) to the equation (31) for n = 3: 

-1 



\wi,W2 = Wi) 



51/(3) 



1 



wi{a[ + al) + a\ 
(1 + 2wlwi] 



V2wibl + h\ 
(1 + 2wlwi)^ 



|0) 



n N 



|0) 



(32) 



where we performed a change of basis in the single-particle Hilbert space, corresponding to 
the following unitary transformation of the bosonic creation operators 

-61 = -!(.; +4) 



h\ = at 



(33) 



V2 



According to the equation (32), when restricted to a invariant subspace under the classical 
dynamics, the SU(3) coherent states are reduced to the SU(2) coherent states with parameter 
72^1. 

Also notice that the state presented in (32) has zero occupation number in the mode 
associated with the operator b\. Therefore, the constraint (30a) is classically equivalent to 
the equation {b^b^) = 0. However, by applying the transformation (33) to the Hamiltonian 
(27), we can easily see that the mean occupation (6363) does not remain zero under the 
quantum evolution of the condensate, considering any state initially unoccupied in this 
mode. Consequently, the subspaces (30) do not have quantum counterparts with identical 
characteristics. However, we can still use the SU(2) coherent states to approximate the 
semiclassical dynamics under these restrictions. This approximation should provide accurate 
results when a similar evolution in the unrestricted space displays irrelevant values of (bibs). 



B. Classical approximation 

In order to establish a criterion for comparison between the semiclassical and quantum 
results, we now introduce a third approach to the bosonic dynamics, which we call classical 
approximation. 

We designate as principal trajectory, indicated by the subindex 'p', the solution of the 
classical equations of motion (6) subject to initial conditions Wp{0) = Wt and Wp{0) = w*. 
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In this case the two vector equations in (6) become redundant, since the solution is such 
that Wp{t) = wl{tf^\ 

The classical approximation to the mean of an arbitrary observable O at the time r is 
defined as follows: 



where \wp{t)) indicates the coherent state parametrized by the principal trajectory Wp(r). 
The classical approximation consists simply in calculating the function On{w*,w) = 
{w\0\w), which represents the normally ordered symbol of the operator O, on the prin- 
cipal trajectory. 

Assuming an initial state Iwi)^^^^^ the classical approximation of (0)(t) is exact in only 
two specific situations when compared with the corresponding quantum results: (i) for 
H G su(n), because in this case |wp(r)) differs from the correct solution of the Schrodinger 
equation by no more than a global phase ; (ii) in the macroscopic limit, given hj N ^ oo"'^. 

Clearly, the semiclassical approximation is more accurate than the classical approach 
(34), since it adds quantum corrections to the classical results. Therefore, the semiclassical 
propagator (5) is also exact for any linear Hamiltonian in the generators of SU(ri) {H G 
su(n)) as well as in the macroscopic limit (A^ — )■ oo). 

Under the restriction H G su(n), every initial condition Wi must provide a trajectory 
with appropriate contribution to the integral (21). Accordingly, the heuristic filter (22) 
must allow the contribution of all trajectories at all instants of time, which it does, because 
\Ksc{w{t),Wi]t)f is constant with respect to t for linear Hamiltonians. 

It follows that the classical and semiclassical approximations to the Hamiltonian (27) are 
exact for x = 0, since in this regime H is linear in the generators of SU(3) (bilinear in the 
creation and annihilation operators). Therefore, the bosonic collisions introduce nonlinear 
terms to the condensate dynamics, whose classical and semiclassical descriptions are not 
complete for a finite number of particles. Consequently, we expect the application of the 
semiclassical propagator (20) to be better behaved for weak nonlinearities (small values of 
x) and large numbers of bosons. 




(34) 
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C. Semiclassical approximation with SU(2) coherent states 



A relevant observable in the condensate dynamics is the population imbalance operator 
Sz, which describes the difference in occupation between the two effectively occupied modes 
in the classical invariant subspace (30a): 

blh - h\h2 



(35) 



Figure 1 compares the semiclassical, quantum and classical evolution of {Sz)/S for 



30, n = —1 and x = ~1) considering as initial state \\/2wi)^ 



SU{2) 



8)su(2) 



. The mean 



of Sz is normalized by the quantity 5* = y so that — 1 < {Sz)/S < 1. For the semiclassical 
approximation we used the SU(2) propagator with 479 initial conditions Wi and limiting 
value A = 10 for the heuristic filter. 
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FIG. 1. Classical (dashed green), quantum (dotted red) and semiclassical (solid blue) evolution of 
the normalized mean of the population imbalance operator Sz for the initial state \^/2wi)g-^^2) ~ 
I tan ^)gjjj2) • The parameters of the Hamiltonian assume the values N = 2S = 30, fl = —1 and 
X = —1- The SU(2) semiclassical propagation was performed with a grid of 479 initial conditions 
and limiting filter A = 10. 



Notice that the oscillations of the classical mean display constant amplitude, unlike the 
semiclassical and quantum results. Although restricted to the SU(2) propagator, the semi- 
classical method shows quantitative agreement with the exact quantum calculations, being 
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fairly superior to the classical approximation, even for a relatively small number of particles. 
In general, the classical and semiclassical approximations are accurate for sufficiently short 
times, but the quality of the semiclassical evolution is obviously higher for longer periods of 
propagation, when the nonlinear terms of the quantum Hamiltonian become important. 




I I I I 1 I 1 I 1 

2 4 6 8 10 

\a\t 

FIG. 2. Semiclassical evolution of {Sz)/S for = 30 (solid blue), = 60 (dotted red) and 
N = 150 (dashed green). The dash-dotted black curve represents the classical approximation, 
which is equivalent to the macroscopic limit N ^ oo. In all results we consider il. = —1, x = ~^ 
and initial state \\/2wi)g^.^, = \ tan f )gu(2)- 



Figure 2 shows the behavior of the semiclassical evolution of {Sz)/S with the variation of 
the total number of particles, forf2 = — l,x = — 1 and initial state |-\/2iy 1)3^^2) = I tan |)gu(2)- 
The results correspond to the SU(2) semiclassical propagator for 30, 60 and 150 particles, 
with A = 10 and about 500 initial conditions in each case. 

Note that the equations of motion (29) and their solutions, including the principal tra- 
jectory Wp{t), are independent of the total number of particles. Therefore, it is easy to show 
that, for a linear operator in the generators of SU(3), the classical mean per particle is also 
independent of A^. Therefore, quantities like {Sz)c/S represent the macroscopic limit of their 
quantum and semiclassical counterparts, since the classical approximation (34) is exact for 

^ 00. 

In accordance with the previous discussion, we included the classical approximation in 
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figure 2 as the macroscopic limit for tlie dynamics of the semiclassical means. Note that 
the semiclassical results quickly converge to the classical curve with increasing A^. Conse- 
quently, we expect the classical approximation to show high accuracy for a few hundred 
condensate bosons, which represents a scenario compatible with usual experiments. How- 
ever, the semiclassical propagators must provide superior results for the mesoscopic dynamics 
when subjected to longer periods of propagation or more intense nonlinear effects. 




-0.2 0.2 0.4 0.6 0.8 1 



FIG. 3. Diagram of contributing trajectories for the SU(2) semiclassical propagator with N = 
30, $7 = —1, X = ~1 Slid initial state l^/lwi)^-^^^^ = \ tan^)^-^^^y Each initial condition Wi 
is represented by a square whose color indicates the period of contribution of the corresponding 
trajectory, according to the heuristic filter (22) for A = 10. 

Figure 3 shows the diagram of contributing trajectories for the SU(2) semiclassical prop- 
agator with N = 30, ri = —1, X = ~1 ^-iid initial state \\/2wi)g-^^^^ = Itanpg^^^)- This 
diagram corresponds to the semiclassical approximation shown in figure 1 and reproduced 
in figure 2. Each square in figure 3 represents an initial condition Wi used in the numerical 
calculation of the integrals (21). The color code indicates the time of contribution of the 
resulting classical trajectories, determined by the heuristic filter (22) with A = 10. 

Notice that the trajectories with the most significant contributions have initial conditions 
Wi centered around w* = tan | ^ 0.29. This initial value defines the principal trajectory, 
whose contribution is among the most important in the reconstruction of the semiclassical 
propagator. Note also that w* is the value that maximizes the Q representation for the 
state |-\/2w*)gyj2)- Therefore, this coherent state is located in the same region of phase space 
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responsible for the most relevant contributions to the initial value representation. 
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FIG. 4. Classical (dashed green), quantum (dotted red) and semiclassical (solid blue) evolution of 
{Sz)/S for N = 30, = —1, X = ~8 and coherent initial state \\/2wi)g^j^^^ = |tan|)gy^2)- The 
semiclassical approximation was performed with the SU(2) propagator, considering a grid of 3781 
initial conditions and A = 18. 

In general, the equations of motion resulting from the Hamiltonian (27) show significant 
changes in behavior for different magnitudes of the ratio x/^ ' ? which represents the 
relative intensity between the quadratic and linear terms of H. The previous examples of 
application of the semiclassical propagator are restricted to small absolute values of x/^j 
since the linear terms are clearly dominant in the dynamics of the condensate. Figure 4 
displays the semiclassical, quantum and classical dynamics of {Sz)/S in a strongly nonlinear 
regime, for = 30, = —1, x = ~8 ^ind initial state \\/2wi)g-^^2) ~ I ^^'^ f )su(2)- 
semiclassical approximation, we employed the SU(2) propagator for a grid of 3781 initial 
conditions and limiting value A = 18. 

Again we see that the amplitude of the classical mean remains constant during the whole 
evolution of the system. Conversely, the semiclassical and quantum results exhibit an al- 
most complete 'collapse' of the oscillations, followed by a partial 'revival' of the amplitude 
value in relation to the classical approximation. Therefore, this example refers to a strongly 
nonlinear and exclusively quantum behavior, described with excellent accuracy by the semi- 
classical propagator. However, note that the number of trajectories required for a proper 
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semiclassical approximation is considerably larger than in the predominantly linear dynamics 
shown in figure 1. As expected, the semiclassical propagator loses computational efficiency 
in nonlinear regimes. 




X y X y X y 

FIG. 5. At the top (bottom) we show the Q representation on the unit sphere related to the 
semiclassical (quantum) evolution of condensate at three different times, for N = 30, = — 1, 
X = — 8 and initial state \^wi) ^^^^^ = \ tan|)gy^2)- 

The phase space corresponding to the SU(2) coherent states may be identified as a spher- 
ical surface"''. It follows that, applying the definition (23) with the coherent states given by 
(32) under the transformation of variables \/2wi = e~*"^tan|, we obtain the Q representa- 
tion for SU(2) in terms of angular spherical coordinates. In this way, we can represent an 
arbitrary quantum state on the unit sphere: 

X = [Q{0, (p) + 1] sin 6 cos (f) 
< y = [Q{e,(j)) + 1] sine sin (f) (36) 
^ ^ = -[Q(^,0) + l]cos^ 

where (j) G [0, 27r) and 6 = [0, it]. Notice that our definition for the variable 6 has its origin 
in the negative z semi-axis. 

In figure 5 we show the comparison between the semiclassical (top) and quantum (bottom) 

18 



Q representations at three different times, for = 30, f2 = — 1, x = —8 and initial state 
|\/2wi)gyj2) ~ l^^^f)su(2)- '^'^^ represented states are in correspondence with the results 
displayed in figure 4. 

At \VL\t = we show the initial coherent state, whose representation is identical in the 
semiclassical and quantum approaches. At the time \Q\t = 3.1, we have the superposition 
of two localized states in phase space ('Schrodinger-cat' state), which is responsible for the 
oscillation collapse in {Sz)/S. At \Q\t = 6.75, we see that the Q function converges again to 
a single location on the sphere. This behavior is associated with the revival of the oscillations 
in figure 4. 

The differences between the quantum and semiclassical representations in figure 5 are 
almost imperceptible, evidencing that the semiclassical approximation accurately describes 
the delocalization and the subsequent relocalization of the state in the phase space. 

D. SU(3) semiclassical propagator 

Although the approximations with the SU(2) semiclassical propagator have shown excel- 
lent accuracy, the SU(3) coherent states are more appropriate to the dynamics determined by 
the Hamiltonian (27). Figure 6 exemplifies the use of the SU(3) semiclassical propagator in 
the evolution of {Sz)/S, for = 30 f2 = —1, x = — 1 and initial coherent state parametrized 
hj Wi = W2 = ^ tan |. In the calculation of the initial value representation we used 35134 
classical trajectories, whose contributions were determined by the heuristic filter (22) with 
A = 10. In comparison with the result for the SU(3) propagator, we reproduce in figure 
6 the corresponding SU(2) approximation and the exact quantum evolution, also shown in 
figure 1. 

As expected, the SU(3) semiclassical propagator is more accurate than the SU(2) ap- 
proximation. The difference between these results comes mainly from the occupation of 
the mode associated with the operator b\. During the considered period of propagation, 
the normalized mean {h\h^)/N grows monotonically until it reaches a value close to 0.04 at 
\VL\t = 6. 

We conclude that most of the inaccuracy attributed to the SU(2) semiclassical propagator 
in figures 1 and 6 is due to the classical constraint (30a), since the SU(3) semiclassical 
approximation is almost exact in the predominantly linear dynamical regime. 
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FIG. 6. Time evolution of {Sz)/S resulting from the SU(2) semiclassical propagator (dashed green), 
the SU(3) semiclassical propagator (solid blue) and the exact quantum calculation (dotted red), for 
A^ = 30, il = — 1, x = ~l and initial coherent state parametrized hy wi = W2 = ^ tan ^. In the 
SU(3) semiclassical approximation we used 35134 classical trajectories, with contributing period 
determined by A = 10. The SU(2) semiclassical curve is the same one shown in figure 1. 



We constructed an initial value representation for the SU(r;,) semiclassical propagator, 
which replaces the search for boundary- valued trajectories by an integral over a set of initial- 
valued trajectories in the doubled phase space. This formulation represents a considerable 
advantage in the calculation of the propagator, since the numerical or analytical resolution 
of a boundary condition problem is typically much more difficult than its initial condition 
counterpart, particularly in systems with many degrees of freedom. Moreover, our method 
allows the factorization of the arrival point wj, as given in equation (20), considerably 
reducing the number of integrations required for a complete representation of the system. 

The semiclassical approach showed excellent accuracy when compared to exact quantum 
results, even for a relatively small number of particles. The efficacy of the semiclassical ap- 
proximation is largely due to the effective heuristic filter, which was able to discriminate the 
trajectories with appropriate contributions to the propagator. The systematic elimination 
of non-contributing trajectories represents a crucial component in the implementation of an 
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initial value representation in the doubled phase space, because it directly determines the 
speed, precision and applicability of the method. 

We tested our semiclassical formula for a triple-well Bose-Einstein condensate in nonlin- 
ear and predominantly linear dynamical regimes. Although the semiclassical propagation 
has been very satisfactory in both situations, the number of initial conditions required for 
an appropriate description of the nonlinear dynamics is significantly higher than in the al- 
most linear case. In general, the computational efficiency of the semiclassical propagator 
is only limited by the required number of contributing classical trajectories. Clearly, this 
number grows with a exponent proportional to {n — 1), the dimension of the subspace w. 
However, we can assume that the required number of initial conditions Wi decreases with the 
total number of particles, since the semiclassical results converge with increasing to the 
classical approximation, which is determined by a single trajectory. Therefore, the SU(?7,) 
semiclassical propagator is a viable alternative in the study of bosonic systems with many 
degrees of freedom and large number of particles, since the computational cost of exact 
quantum methods typically grows as a polynomial in of order proportional to n. 

Finally, we would like to point out that the formulas (20) and (21), the main results 
of this paper, can be easily extended to other classes of coherent states, such as the usual 
harmonic-oscillator coherent states. Thus, this work also represents an alternative to previ- 
ously published semiclassical methods. 
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